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SUMMARY 


The role of fluid-forces acting on the blades of an axial turbo-rotor with 
regards to whirling is analyzed. The dynamic equations are formulated for the 
coning mode of an overhung rotor. The exciting forces due to the motion are 
defined through a set of "rotor stability derivatives" (R.S.D.'s), and analyti- 
cal expressions of the aerodynamic contributions are found for the case of 
small mean stream deflection, high-solidity and equivalent flat plate cascade. 
(Torque-whirl and tip-clearance effects can also be included in the R.S.D.s). 

For a typical case, only backward whirl is indicated if the phase-shifting of 
the rotor wake effect is ignored. A parametric study of the dynamic stability 
boundary reveals that a reduction in blade stagger angle, mass-flow rate, fluid 
density and an increase in stiffness and external damping are all inducive for 
improved stability. An optimum overhang distance of the rotor from the bearing- 
support can also be found. Finally, when two or more opposing whirling mecha- 
nisms are present, mutual annihilation is possible by making a certain whirling 
group "b 0 " very small. This concept can be useful in the preliminary design 
stage or for later improvements. 


INTRODUCTION 


An increasing number of severe failures in high speed compressors and 
turbines in recent years are attributable to whirling instability. Ehrich 
(ref. 1), Shapiro and Colsher (ref. 2) and many others have reviewed the de- 
stabilizing mechanisms. To our knowledge, only a few studies have been con- 
ducted regarding the fluid-dynamic forces acting on the blades themselves. 
Alford (ref. 3) dealt only with the periphery of the rotor, namely, the tip- 
clearance and labyrinth seal forces. Reference 4 includes various other mecha- 
nisms acting between the rotor-circumference and its casing. Bousso (ref. 5) 
and later Vance (ref. 6) analyzed the consequences of the tilting of load- 
torque, but did not account for other effects produced by the same fluid 
forces. Trent and Lull (ref. 7) noted the important analogy between turbo- 
rotor whirl and propeller whirl. Ehrich (ref. 8) considered the fluid react- 
ions based on Euler's turbine equations but his dynamic equation was later 
found to be in error. 

In this paper, we adopt the usual modelling of an overhung rotor in its 
'coning mode' to set up the dynamic equations of motion. The forcing terms are 
from the self- induced aerodynamic forces due to the motion, and the whirling is 
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thus reduced to an eigenvalue problem much like flutter. The analysis follows 
the basic idea of propeller whirl but the details are carefully reworked for 
the cascade configuration in turbo-rotor systems. Our emphasis is to lay down 
a proper framework and outline the methodology, so that only systematic 
refinements and extensions are needed to yield results of practical 
significance . 


SYMBOLS 


a^,b^ displacements of shaft-center in 
y- and z-directions respectively 

diagonal damping matrix with 

elements c.e^ 
d 

C complex wake function 

coefficient of lift 

C Rotor Stability Derivative 

a B (R.S.D.) 

c chord length 

2 

D = 2?K Q (e/r 2 ’) 

e equivalent overhang shaft length 

F real part of wake function C 

G imaginary part of wake function 

C 

H = J/l, ratio of inertias 

I diagonal inertia matrix with 

elements I 
2 

I = Mr 2 , moment of inertia of rotor 

about pivotal axis 

i,j unit vectors in axial flow 

direction and opposite to 
rotational speed, respectively 

J polar moment of inertia of rotor 

K diagonal stiffness matrix with 

— elements K 

K,k 2 equivalent torsional spring 

k cascade interference factor 

o 

K = w R/V , reduced structural 

9 -9 a 

frequency 

k reduced aerodynamic frequency 

Z elemental lift 


M mass of rotor 

2 2 

M' = ttR pV(l - n ), annular mass flow 
. o 

rate 

m elemental moment 

r radial location 

2 2 1/2 

r 2 (r 2 + e ) , radius of gyration 

of rotor about pivotal axis 

R rotor tip radius 

s gap between adjacent chords 

T kinetic energy 

t time 

V potential energy 

V absolute velocity 

W relative velocity 

3 angle of velocity component with 

respect to axial direction 

y blade stagger angle 

A change 

5 angle of attack 

n = r/R 

n hub-to-tip radius ratio 

o 

0 amplitude of pitch angle 

o 

0^ coning angle 

Q. azimuthal location of i th blade 
from horizontal axis 
5 2 

k = ttR p(l - n )/l, ratio of moment 

of inertia°of annular cylinder of 
air to that of rotor 

p fluid density 

a = c/s, solidity ratio 
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4 > = 

V/QR, tip mass-flow coefficient 

comp 

compressible 

^o 

amplitude of yaw angle 

ef f 

effective 

CO 

angular velocity vector 

inc 

incompressible 

O3 0 = 

/k7i 

Im 

imaginary part 

Subscripts, prefixes, etc. : 

Re 

real part 

o 

steady value 

s 

inertial reference frame 

1,2 

inlet and exit stations 

x,y,z 

in the direction of x,y,z axes 

1,2 

forward and backward whirl 

w 

wake contribution 

1,2 

refers to (Q + to) in context 
with reduced frequencies for 

0 

tangential direction 


whirl reference value 

o 

average value 

a 

axial direction 

(’) 

perturbed value 

( ) 3 

aerodynamic contribution 

( ) 

matrix notation 


THE ROTOR MODEL 


A rotor on a flexible shaft is equivalent to a gyroscope tied to a rotat- 
ing spring. Figure 1 shows the real rotor shaft system and the idealized model 
considered in this paper. This simplest model focuses on the 'coning mode' and 
has all the characteristics necessary to produce the phenomenon under study. 
Some of the implications and further assumptions are: 

1. The shaft elastica remains in one plane due to infinite stiffness in 
torsion. 

2. The angular tilt of the rotor axis essential for the coning mode is 
coupled to the radial deflection of the center of the rotor. 

3. The cantilevered shaft is approximated by an equivalent, torsional spring, 

and an equivalent, rigid, massless shaft of length, lj_. The notation 
K = ki and e = 1]_ is used in the main text. (See reference 9 for 
inadequacies in this model.) 

4. The bearings are rigid and frictionless. 

5. The center of mass of the rigid rotor is at its geometric center, i.e. 
there is no imbalance or eccentricity. 

LAGRANGE'S EQUATIONS OF MOTION 

The motion and instantaneous position of the rotor can be completely 
specified by the Euler angles ip , 0 and Qt which are, respectively, the yaw, 
pitch and azimuth angle as shown in figure 2. For small angular deflections 
i > , 9 and constant rotational speed we have 

• • • 

a) = n + ipd , 03 - -0 and to = ip (1) 

x y z 

The dynamic equations are derived by using Lagrange’s equations: 
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d_ 

dt 


3L 


( 2 ) 



9q i 


9D 




where, 

It It 

L = Lagrange function = T - V = — u — q • K • 

1 m 2 — 2 — — 

D = Dissipation function = — q • C • q 
= Generalized force = 6W/ 6q^ 

6W = Virtual work 

and q^ = Generalized degree of freedom (d.o.f.); (q^ = ip, q^ = 
After simplifications, the results may be written as: 



q 


e) 


(3) 


In particular, note that I is the moment of inertia of the rotor about the 
"pivotal axis" and c^ is the "equivalent viscous damping coefficient" to ac- 
count for the damping effects assumed to be representable by an effective force 
applied at the tip of the shaft. 

According to the principle of virtual work, 5W = M z + M 60. It fol- 
lows that = M z and = My, i.e. the generalized forces are the pivotal 
moments M z and My produced due to the fluid-dynamic forces acting on the rotor. 
These moments are found in the next section to be functions of both the gener- 
alized d.o.f., q^, and their derivatives. Thus, besides the anti-symmetric 
gyroscopic coupling due to + JQ there also exists coupling due to the fluid- 
dynamic forces. Such external couplings, if anti-symmetric in nature, as is 
well known, can produce unstable response. 


THE ROTOR STABILITY DERIVATIVES (R.S.D.'s) 


The generalized forces caused by small perturbations of the d.o.f. will 
be expressed as a set of 'rotor stability derivatives', abbreviated as R.S.D.'s, 
in this paper. For a systematic treatment of the aerodynamic terms, the anal- 
ogy of the present problem to propeller whirl treated in References 10 and 11 
is obvious. However, we have to make judicial adjustments for the cascade 
configurations. As a first approximation, we consider here only the quasi- 
steady effect of inviscid, potential fluid flow on an axial turbo-compressor. 


The quasi-steady lift for an airfoil in cascade depends on the instantan- 
eous angle of attack and is given by 


1 - \ 


C£ eff 


c Ar 


(4) 
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In order to find the effective reference velocity, W e ff, and the coefficient 
of lift Cj>, e f£, consider the displaced rotor as shown in figure 3a. A circum- 
ferential cut AA at any radial location unfolds a cascade of airfoils which 
can be represented simply by equivalent flat plates as shown in figure 3b (see 
ref. 12). It is well known that, for a cascade of high solidity, i.e. small 
gap-to-chord ratio (s/c $ 0.7), the angle of the relative outlet velocity with 
respect to the blade is independent of the relative inlet flow angle and is 
equal to the equivalent flat plate stagger angle. Hence, assuming that the 
rotation of the blade due to yaw and pitch is very small as compared to^its 
stagger angle, y, the perturbed and steady relative outlet velocities, W 2 ' and 
^ 2 > are approximately parallel to each other. Next, assuming constant radius 
cylindrical stream surfaces, the equation of continuity gives the equivalence 
of the axial components of the absolute outlet velocity in the steady and per- 
turbed states, i.e. = ^ al ~ V a . The velocity triangles in figure 3c are 
shown using these assumptions. The velocity triangles are crucial to the 
development of the R.S.D.'s. (See ref. 9 for further discussions.) 

Defining the reference velocity in cascades by 

W = \ (W. + W ) (5) 

00 2 1 2 

we get + ^ , 1 

AW = W* - W = vi + — (u + vtan 3„)j (6) 

OO CO CO 2 2 

where, u and v are the perturbation velocities of the blade in the plane of 
the rotor disc and normal to it. From kinematical relationships, we can easily 
write 


u - — eip sin fit + e 0 cos fit 

v = (i^fi - 0) r sin fit - (0fi + ip) r cos fit 

The angular tilt of the rotor axis and the component of AW^ normal to W a 
combine to give an effective change in the angle of attack, 

A<5 = 6 + 6 2 


(7) 


( 8 ) 


where, 

6 1 = ^ 

sin fit - 0 

cos fit 



< 5 „ = - 

v (fir + V 

tan - V ) /2W 2 

+ i 


2 

a 

2 01 00 

2 

Using the coefficient 

of lift for 

■ a flat-plate cascade ( 


c„ 

= 2n k sin 

5 



£ cc 
ef f 

0 

eff 


with 

<5 ^ 

= 6 + AS 




eff 

0 



and putting 

W 

eff 

= W + AW 

CO 








where 

AW 

= component 

of AW OT parallel to 

Wco , 


( 9 ) 



we get the perturbed component of lift from equation (4) as follows: 

l' 2 

— - — r— = W cos 6 A6 + 2 W sin 6 AW 

irpk c Ar oo o °° o 

o 

2 2 

+ 2W cos 5 AW A6 + sin 5 AW + AW A5 . (10) 

oo o o 

The first term provides the leading contribution and for lifting cascades (i.e. 
6 Q ^ 0) the second term gives an additional small effect. Equation (10) is a 
general expression for the perturbed lift component; it can be written for the 
i th blade by replacing fit by fit + 8^, then summed over all the blades and 
finally integrated from the hub to the tip radius giving the total perturbed 
lift. 


We now specialize equation (10) for a non-lifting cascade, i.e. <$ 0 = 0 and 
uniform axial inlet field with no swirl, i.e. = 0. It turns out that even 

in such a simple case significant forces may be produced. 


Using the simple relation for the cascade interference factor, k Q 

k = ' s _1_ 

o TT c COS Y 


(ID 


for a high solidity cascade with (s/c) 0.7, (see ref. 12), the following 

result for the non-dimensional force components can be obtained: 


where 


= C ip + C 
MV y y 

a ip J 


\p R 
V 

q \ a 
0 R 


— r^— = C 0 4- C 
MV z z \ V 

a 9 r \ a 


+ C 


+ C 


c - c = 1 + (l + n 2 )/(4(j> 2 ) 
z e 

C = C = - e/2R 
y z 
q r 



c = -c = (i + n 2 )/(44> 2 ) . 

y z o 

r q 


Similarly, integrating the elementary moment components 


m = - r l tan 3 

y y co 

and 

m = r i tan g 
z z 00 


we get the non-dimensional moment components as 


( 12 ) 


(13) 


(14) 
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where 


mVr ' c m * + c m (1^) + C m (f^ 

a <p q \ a / r\a 


M 


M’V R n Q 

a 9 


= C 0 + C 


„ (P) + (F) 

r ' a / q\a/ 


( 15 ) 


r ry « / q 

C = -C = - (1 + n )/2(j> + (1 + t ) + n >/6(t> 

"V n 0 L ° o o 

c = -C = ff (i + n 2 )/<f> 

m n 4 R o 

q r 

C = C = -(l + n 2 + nS/6<j> 2 . 

m n o o 

r q 


(16) 


Analogous to the propeller-in-yaw studies, for example, by Ribner (ref. 
perturbed axial force and torque about the rotor axis are zero. 


13) the 


The coefficients C, 


etc. 


, \jy , , Ojjj , vjjjj , , ci_l. as defined by equa- 

tions (12-16) and which denote increases in fo?ces/moments due to unit increases 
in the d.o.f. or their derivatives are termed as the 'Rotor Stability Deriva- 
tives' in this paper. Note that in this formulation the R.S.D.'s are independ- 
ent of the number of blades and solidity ratio and remind one of actuator disc 
type results. In practical applications reliable determination of the R.S.D.'s 
is, of course, imperative. 


APPROXIMATE CORRECTIONS FOR COMPRESSIBILITY AND WAKE EFFECTS 


For compressible flows, a first-order correction on the pressure distribu- 
tion can be obtained by using the expression given in reference 12, page 61. 

For the case of small lift and low subsonic Mach numbers it reduces to the well- 
known Prandtl-Glauert Rule: 


(C 0 ) = (c 0 ). / 

£ comp l me T 


M 


where M is the local Mach number at any radial location n. The concept of an 
effective Mach number, M e ff was used by Ribner (ref. 13) as an overall correct- 
ion to the R.S.D. So long as M is small, we find only small changes in the 
R.S.D.'s and the stability boundaries are not significantly affected. 


The unsteady rotor wake can be taken into account by the use of complex 
rotor wake functions; see Houbolt and Reed (ref. 10) in the propeller whirl 
problem. Instead of the well known Theodorsen function for an isolated airfoil, 
one needs the oscillating cascade wake functions, C, for example, by Whitehead 
Cref. 14). In cascades the wake function is a parameter of not only the 
reduced frequency, k, but also of solidity, blade stagger angle and the phase 
of vibration between adjacent blades, $. Hence, 
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and 


F = Re { C } = fn (k,s/c, y , $) 
G = Im { C} = fn (k,s/c, y , $) 


( 17 ) 


Analysis done parallel to reference 10 leads to the following important 
conclusions for multi-bladed cascaded rotors (see ref. 9): 

1. Two reduced frequencies are involved in whirling: 

k- = (£2 + id) c/2W (18) 

1 , Z — °° 

2. If the whirling frequency ui is much smaller than the rotor 
frequency Q, then the effect of wake is to turn the resultant 
force and moment clockwise by approximately: 

S* = tan _1 (G/F) (19) 

where 

F-| j F(k x ) + F(k 2 ) j 

G = j j G(k x ) + G(k 2 ) \ (20) 

without much change in the magnitude of the vectors. 

3. Two new R.S.D.'s solely attributable to this effect are added, 


namely, 

c w = 

-c w 

= c 

tan 6* 

O 

CD 

II 

Z *l> 

c w 

= c 

tan 5* . 

m 9 


m. 



Hence, the contribution to the effective tangential force due to the angu- 
lar tilt of the rotor, which is important for whirl, appears from two sources: 

i) C : the quasi-steady moment R.S.D. 
m , 

ii) C w : the unsteady rotor wake force R.S.D. 

Z | 

Ehrich's analysis (ref. 8) did not include either of the above two mechanisms. 

The above analysis which was for a compressor or a negative work turbo- 
machine can be repeated for a turbine or a positive work turbo-machine. Though 
the previously derived R.S.D.'s are not applicable for thick, curved turbine 
blades with high flow deflection, it can be shown that the direction of the 
net force and moment does not change. 

For lifting rotors, two other important mechanisms discussed in the liter- 
ature are Alford's tip-clearance effect (ref. 3) and Bousso's torque whirl 
effect (ref. 5). Both can be represented by an equivalent tangential force due 
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to a radial displacement and contribute thus only to C (= -C ). In ref. 9 

Z V 

the conversion to R.S.D. is explicitly given. ip J Q 


DYNAMIC STABILITY ANALYSIS 


by 


The generalized forces are related to the fluid-dynamic forces and moments 


Q = M = eP + M 

1 z y z 

Q„ = M = eP + M a 

2 y z y 


(22) 


It is advantageous to write the equation of motion (eqn. (3)) in non- 
dimensional form. Defining a new independent variable 


T = V t/R 
a 

and denoting derivatives with respect to r as superscripted primes, 

.2 


1 

0 


°U*" 

0" 


+ 


I'D H/<bW0 ’ 
)~H/ <p DH0 ' 


t- 


(23) 


(24) 


where the non-dimensional forces are given by 


with 


f =a 0 + a 0' + b 0 + b 0' 

ip o 1 o 1 

f - a 0 + a 0' - b 0 - b 0' 

0 o 1 o 1 


a = ^ C + C 
° R ^ m 0 

a = f C + C 

1 R Yq m r 


b = - C 
o R y c 


- C 


m , 


b l " ^ C 


„ _ - C 

R Vr m n 


(25) 


(26) 


Hence, the homogeneous equation is 


i o (V-i 


i D-Ka, 


H/ 0-icb 


:o i)M \ -H/0+tcb^ D-<a^ j 


ir< 

0’ 1 


( K A" Ka 

) 9 o 
I tc b 


~ Kb o 1 
K g-<a i 

0 O' 


f'f'l 

le 


0 

0 


(27) 


The anti-symmetric off-diagonal terms (H/0 - Kb^) and (-<b Q ) are important 
to the self-excited whirl phenomenon as noted earlier. 
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In order to find the conditions for which the system becomes neutrally 
stable, sinusoidal motion for \p and 8 is specified: 


* = K 


iwt 

e 


0 


„ iwt 
6 e 
o 


( 28 ) 


where w, the whirling frequency is assumed positive. 


To avoid a frontal attack of the eigenvalue problem with a large number of 
parameters, the following approach has proved to be more fruitful. Let D in 
equation (27) be replaced by D which denotes the damping required for neutral 
stability. Then the substitution of equation (28) in equation (27) gives rise 
to two algebraic equations in terms of the two unknowns w and D. Equating the 
real and imaginary parts to zero gives 


and 


A = 


2K„ 


b + \fo - 4d J 


£ = 


D 


2K 0 (e/r’) 2 2K 0 (e/r’) 2 


(H/<J)-Kb )Kb Q 

<a + — 2 2 

' <(l-X Z )-Ka 
y o 


(29) 


(30) 


where A = w/w , non-dimensional whirling frequency 

0 

t, = c^/2Mwq, fraction of critical damping required for 
neutral stability 

b = (H/<f>) 2 + 2K 2 - 2 k (a + Hb /<j>) + K 2 ( a 2 + b 2 ) - 2DKa, + D 2 
0 o I 11 1 

d = (K 2 - k 3 ) 2 + (Kb ) 2 . (31) 

0 o o 

2 

It is necessary for neutral stability that A >_ 0. Further, only the positive 
root for A is taken and (b^ - 4d) > 0 should be satisfied. 


The equations (29) and (30) are non-linear in A and 
^-values of interest, the second order appearance of £ or D 
ly neglected. Letting b* denote the value of b with D = 0, 
mate of A can be obtained, namely, 

1 


but for the k- and 
in A^ can be usual- 
a very good esti- 


A 


1,2 


2K 


b* + 




4d 


(32) 


9 


Substitution of these A^ , A 2 into equation (30) gives the corresponding or 
7 ^ 2 • If the actual damping ^ is larger than the damping required for neutral 
stability, £, we conclude the system must be stable. Otherwise, the system 
is unstable and the amplitude increases unless limited by non-linearities. 
Examination of X (eqn. (30)) reveals that since Ka-^, occurring due to aero- 
dynamic damping, is always negative, a 'conservative estimate' of £ can be 
obtained by dropping Ka^ term. Further analysis shows that generally the term 
(1 - A“) changes sign for the two roots of equation (32). A conservative esti- 
mate of the needed £ is thus proportional to b Q , the whirling group. Hence, 
the larger the magnitude of b Q the larger is the value of damping required to 
prevent whirling. 
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NUMERICAL RESULTS FOR AN EXAMPLE AND DISCUSSION 


Table I gives typical values of design parameters used for numerical anal- 
ysis. They are believed to be representative of small, high-speed turbo-rotors. 
The high value of the hub-to-tip radius ratio is tailored to approximate our 
assumption of constant radius cylindrical stream-surfaces. 

The R.S.D.'s are shown in figure 4 as functions of the tip-mass flow coef- 
ficient, <f > . Since <j> = cos y , we notice that R.S.D.'s increase with the blade 
stagger angle. The whirling R.S.D.'s C and Cy (which is proportional to 
Cy.) are relatively larger as compared to the remaining R.S.D.'s. Note that 
C also depends on tan 5*, where 6* is the wake angle given by equation (19). 
The numerous simplified assumptions under which these R.S.D.'s were formulated 
constrain them to only a small range of blade stagger angles. 

It is instructive to examine the effect of only the gyroscopic coupling 
i JT2. Thus, neglecting all the external aerodynamic forces and the viscous 
damping, a plot of X versus the non-dimens ion less rotating speed, f2 / tog » as giv- 
en by equation (29) is shown in figure 5. At any given rotor speed, the two 
whirling frequencies X-j_ and ^ , corresponding to forward and backward whirl, 
respectively, drift further apart from each other as the parameter H, the ratio 
of moments of inertias of the rotor is increased. Also, at higher rotor- 
speeds the backward whirling frequency is very small as compared to the forward 
one or the rotor speed itself. When a slight damping is present both these 
modes die out. However, without the f luid-dynamic forces, the gyroscopic 
action cannot in itself lead to a divergent type whirl instability because the 
net energy input to the system is zero. 

The inclusion of the R.S.D.'s through k in the whirling frequency expres- 
sion has an imperceptible effect on X itself. Hence, the forced whirling fre- 
quencies are approximately the same as the natural whirling frequencies. Figure 
6 shows X vs. 1/Kq with <f> as parameter. Kg is the reduced structural frequency 
and is a convenient independent variable. A tendency to converge to the 
natural frequency to g is observed as the mass-flow coefficient is increased. 

The stability boundary plots ? vs. 1/Kg are shown in figure 7 for various 
parameters. In figure 7a where <f> is a parameter and wake effects are not in- 
cluded the damping required for forward whirl is negative, i.e. < 0. Hence, 
the quasi-steady analysis predicts no forward whirl. On the other hand, back- 
ward whirl can occur and rotors with small cf> and large 1/Kg values are more 
susceptible to whirl. Other studies in ref. 9 indicate that increasing the 
fluid density or Mach number should decrease stability in general; whereas in- 
creasing the rotor-to-bearing distance has a stabilizing effect as in figure 
7c. This is attributable to the increase in aerodynamic damping term, <ap and 
also due to the form of the fractional damping ? assumed. The effect of un- 
steady rotor wake in terms of the wake angle <5* is shown in figure 7b. In iso- 
lated airfoils the calculation of 6*, say, from Theodorsen's function always 
gives negative 5* for all values of the reduced frequency 0 ^ k < °°, and, hence, 
has a stabilizing effect (ref. 10). In cascades, the cascade Theodorsen func- 
tion C may give positive values of 6*; moreover, since C is itself a function 
of the whirling frequency a more involved iterative procedure is needed to 
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solve for X and 4 . In figure 7 b only arbitrary values of <5* are used to bring 
out the effect of the sign of 5*. In the actual calculation of 6* if << fi 
is indicated, then k = £2c/2Wco can be used for F and G. 

Finally, the importance of other competing mechanisms can be assessed. 

For example, the equivalent tangential force contributions to represent the 
effects of tip-clearance and load-input torque misalignment were added to Cy„ . 
Both these mechanisms being absent for a non-lifting cascade, in reference 9 
we have recalculated the R.S.D.'s and the torque coefficient for an assumed 
case of 6 0 = 5°. The resultant stability plot is shown in figure 7d with e/R 
as parameter. This plot demonstrates the possibility of forward whirl for 
certain values of e/R and also shows that backward whirl can be suppressed. 
According to refs. (3) and (5) both these load-dependent mechanisms give rise 
to forward whirl for aft-pivoted compressors. In our example, unless modified 
by the wake effect, only backward whirl is possible. The unsteady rotor wake 
effect depends on the sign of 6* and no conclusion is possible pending addi- 
tional investigation. The direction of whirling which is observed in a practi- 
cal turbo-rotor is the resultant of all the mechanisms. Since Z “ I b Q | , it 
seems evident that an astute aerodynamic design should aim at making b Q small, 
thus leading to the total elimination of any undesirable whirlT 

CONCLUSIONS 


1. The propeller-whirl analogy can be directly applied to the whirl of axial 
turbo-machinery rotors with modifications for the cascade configuration. The 
basic driving mechanism of orthogonal moments is seen to remain the same. 

2. A set of rotor stability derivatives (R.S.D.'s) for turbo-machinery rotors 
has been suggested. They are fundamental to the dynamic behavior. Express- 
ions specifically derived for the coning mode of an aft-pivoted compressor 
with equivalent flat-plate, high solidity cascade are given. For practical 
applications, experimental verification of the expressions would be 
invaluable. 

3. Quick survey of the dynamic stability as affected by various parameters is 
made possible by introducing artificial damping for neutral stability. 

4. Considerations based on the R.S.D.'s suggest only backward whirl of compres- 
sors and turbines. The role of unsteady rotor wake is outlined but the net 
effect awaits more analysis. 

5. The dynamic stability boundary plots reveal that stability can be increased 
by reducing the blade stagger, mass flow rate, density and increasing the 
stiffness of the shaft and external damping. Though an increase in over- 
hang-to-tip radius ratio increases stability, the effective stiffness is 
decreased. Hence, an optimum e/R can be found. 

6. Various mechanisms can be mutually annihilated by making the whirling group 
b Q as small as possible. This goal may be achievable at the preliminary 
design stage or later by modification of the rotor geometry, flow 
distribution or the structural properties. 
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TABLE I.- TYPICAL VALUES USED FOR NUMERICAL ANALYSIS 


Mass of the rotor M = 21.89 kg (1.5 slugs) 

Radius of tip of the rotor R = 0.305 m (1 foot) 

Hub-to-tip radius ratio . 0 = 0.75 

r o 

Polar radius of gyration-to-tip radius ratio of the rotor r^/R = 0.6 

Radius of gyration of the diametrical moment of inertia- 

to-tip radius ratio t^/R = 

3 3 

Standard density of air P Q = ^'^5 kg/m (0.002378 slugs/ft ) 

* . 

Shaft overhang-to-tip radius ratio e/R = 0.3 


In figures 7c and 7d, e/R is a parameter. 



Figure 1. The rotor model. 
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Figure 5. 


The natural whirling frequencies. 



Figure 6. The forced whirling 
frequencies with 
mass- flow coefficient 
as parameter. 



Figure 7a. Stability boundaries with 
mass-flow coefficient as 
parameter. 





Figure 7b. 


Effect of wake-angle on Figure 7c. Effect of shaft-overhang 
the stability boundary. length, on the stability 

boundary. 



Figure 7d. Inclusion of other effects on the stability boundary 
with overhang shaft length as parameter. 


